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q ' Abstract 

■ In many radar scenarios, the radar target or the medium is assumed to possess randomly 

varying parts. The properties of a target are described by a random process known as the 
spreading function. Its second order statistics under the WSSUS assumption are given by the 
scattering function. Recent developments in the operator identification theory suggest a channel 
sounding procedure that allows to determine the spreading function given complete statistical 
knowledge of the operator echo. We show that in a continuous model it is indeed theoretically 
possible to identify a scattering function of an overspread target given full statistics of a received 
echo from a single sounding by a custom weighted delta train. Our results apply whenever the 
O ■ scattering function is supported on a set of area less than one. Absent such complete statistics, 

we construct and analyze an estimator that can be used as a replacement of the averaged 
, periodogram estimator in case of poor geometry of the support set of the scattering function. 

> 

-rf '. 1 Introduction 

CO 

In the classical delay-Doppler radar system, the echo y(t) that is reflected from a target can be 
\Q • expressed as a superposition of time- frequency shifts of the transmitted waveform x(t), that is, 

o 
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X 



y(t) = J J T)(t, 7 )M 7 T T x(t) dr d 7 , (1) 

where T T x(t) — x(t — r) is a time shift operator, M 7 x(t) = e 2mlt x(t) is a frequency shift operator 
and r\ is the spreading function of the target. If the target has randomly varying components, then 
t](t, 7) is assumed to be a random process. It is often assumed that the radar environment satisfies 
the wide-sense stationarity with uncorrected scattering (WSSUS) assumption [Bel63l IVTOlj . In 
this case, rj(r, 7) satisfies the relation 

E { V (t, T )>(t', 7')} = S(t - t') 5(7 - 7') C(r, 7), (2) 

where C(i, 7) is called the scattering function of the target. 

A problem in radar is to determine C(t, 7) from the echo y Grc68 ( Har02 . A standard approach is 
to view the scattering function C(t, 7) as power spectral density of the 2-D stationary random process 
known as the time-variant transfer function. This gives rise to the averaged periodogram estimator 
of C(t,j). One of the significant restrictions of this approach is the requirement to assume an 
underspread target, that is, the support of C(t, 7) must be contained in a rectangle of area less than 
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one. Building up on the recent results in operator identification |P W061 IKP05I IPRT08I Ffa08l IGP08j . 
we show that it is the area of the support of C(t, 7) itself that matters, and not the area of the 
bounding box. This fact calls for reclassification of the channels that would previously be considered 
overspread due to the non-rectangular geometry of the set suppC(i,7). 
In Section [3] we prove the following theorem. 

Theorem 1. Let C(t, 7) be the scattering function of a radar target. Let the support of C(t,"f) 
have area M < 1. There exists T > 0, a natural number J and a fixed J -periodic sequence Ck, 
such that we can identify the scattering function C{t,^f) from the received echo y(t) = Hx(t) where 
x(t) = J2kez Ck ^ kT (^) * s the sounding signal. 

In the tradition of |PWnt) . we provide an explicit formula (fT2"|) for the reconstruction of the 
scattering function that is supported on a possibly overspread domain of the time-frequency plane 
using a specially constructed delta train as the sounding signal. A procedure for the estimation of 
the scattering function is naturally derived from this reconstruction formula. 

The paper is organized as follows. Section [2] gives an overview of classical delay-Doppler radar. 
In Section [3J we describe the target identification problem and prove the main result, ITheorem II 
We suggest and analyze an estimator in Section [H followed by our conclusions in Section [5] 



2 Overview of radar 

The classical scenario in a delay-Doppler radar system is that a testing signal x(t) is transmitted. The 
signal x(t) might be a short pulse, as well as a wideband linear chirp, coded waveform, pseudonoisc 
sequence, etc. [SkoSO, TayOl]. 

The echo from a point target at distance d and traveling at constant speed v, is considered to be 
of the form 

y = m x{t-t a )e 2 ^ 1 

where to — 2d/c, c is the speed of light, 70 is the Doppler shift , and 770 is the reflection coefficient 
that depends on the distance and the speed of the point target. 

As is customary and indicated in ([TJ , we model the echo realized in a complete environment as a 
continuous superposition of time- frequency shifts of the transmitted signal x(t) |Gre68j . Equivalcntly 
to |T]), from the time-varying impulse response 

h(t,r) = J n{ T , 1 )e 2 ^< t d 1 (3) 

we obtain the linear time-varying representation of the channel 

y (t) = I h(t,r)x(t - r)dr. 



Certain characteristic features of a target or a radar scene, for example, sea clutter [Gaa68, VT01| 
IDT07 . FBMC05], are often modelled as random processes. In this model, h and r\ are assumed to 
be a pair of random processes related via the one-sided Fourier transform ([3]) . 

Let the quantity 

R h (t, S ;Z, V )=E{h(t,0*h(s, V )} 

denote the autocorrelation function of h, where E{} denotes expected value. The usual WSSUS 
assumption in radar applications is that h(t, r) is a) zero-mean in both variables, b) wide sense 
stationary (WSS) in t, and c) has uncorrelated scatterers r, which reads 

Rh(t,t + T;£,ri)=P h (T,t)8(t-rj). 
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Then the scattering function C(t, 7) defined in ^ satisfies 

C(r )7 ) = J P h (t,r) e 2 ^ dt. 

C(t, 7) is the previously mentioned scattering function associated with the channel, and Ph(t,r) 
is the channel's autocorrelation function, or ACF. The computations above and in the remainder 
of this paper involve delta distributions, thus the convergence and equality of the integrals must 
be understood in a weak sense. We preserve current notation for clarity and relegate a rigorous 
treatment of this matter to [PZnt] . 

In radar applications, it is usually assumed that the scattering function C(t, 7) of a target has 
its support on a time-frequency rectangle [0,0] x [—51/2,51/2]. The degree of dispersion of the 
echo is quantified by the spreading factor 510. A target is said to be under spread if 510 < 1 and 
overspread if 510 > 1. In recent years it has become increasingly clear that the correct criterion is 
that the area of the support of the scattering function must be less than one [DSBSIO] . We show 
that it is indeed correct, and require only that supp C(t, 7) is compactly supported within some 
rectangle [0,0] x [—51/2,51/2] without demanding 510 < 1. A classical identification problem in 
radar and communications is to estimate the scattering function of a given randomly varying target, 
or, equivalently, an linear time variant random channel. The general approach is to transmit a signal 
x(t), once or multiple times, and to construct an estimator C(t, 7) for the scattering function C(t, 7) 
using the received echoes. Obtaining measurements in this way is referred to as channel sounding 
in communications. 

3 Scattering function identification 

Let rj be the spreading function of a radar target such that supp C(t, 7) C [0, 0] x [—51/2, 51/2], with 
510 possibly larger than 1. Let the set suppC(i,7) be covered by the translations of the prototype 
rectangle with some J prime number, as follows: 

J 

suppC(i, 7 ) C (J R+ ( aj T,bjB) 

where aj,bj £ Z index the cover, T — 0/51, B = 51/ J0. This implies JBT = 1. Note that other 
bounding boxes, for example, a first quadrant box [0, 0] x [0, 51], can be accommodated by a simple 
translation; the case JBT < 1 can be padded with empty boxes, if necessary. The theory of Jordan 
domain shows that for any compact set M of area less than one the prime number J and the covering 
bj}j =1 can always be found. Given 77, we define patches of -q to be 

Vj(t,7) = XR(t,l)v(t + ajT,1 + bjB), 

and their Fourier transforms 

h j (x,t)= J Vj (t n )e 2 ^ x d 7 . 

We can now reconstruct 

J 

v(t, 7) = ~ a J T ' t ~ 
3=1 
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and 



h(x, t) = f Vj(t - ajT, 7 - bjB) e 2m ^ x d 7 
= STe^ B h J (x,t-a 1 T). 



3 = 1 

For each t fixed, rjj(,t) is supported on [0,5). Therefore, for any r we have the orthonormal 
expansion 



(4) 



In particular, setting t = Tj(i) = t + ajT, we obtain 

,/ 

7) = XI - a j T , 7 - bjB) 

J 



j) 

For the random r/j 's, we have 

E {»&(*, 7)*»?/(*, A)} 

= J (t - s + (oj - h)T) 5{ 1 -\+ (b 3 - le )B) x«(t, 7)Xfl(*> A)C(t + 0j T, 7 + 
= ,5 (t - s + {a 3 - t t )T) S (7 - A + (6, - 7/ )B) Xfl (t, 7)Xfl(* + (<*j - f/)T\ 7 + (&j - 7«)B) 
x C{t + ajT^ + b J B). 

But this is zero unless Oj = tt and 6j = 7^, that is j — £, in which case we have 

E 7)*^(«, A)} = 7,<K* - *)«(7 - A)C (t + 0j -T, 7 + • 
In particular, defining Cj(t, 7) = Xfl(^ 7/C(* + a j^\ 7 + bjB), we have the expansion 

J 

C(t, 1 ) = Y,C j (t-a j T, 1 -b j B). (5) 

Similarly, for the autocorrelation of h's we have whenever j ^ I 

E{h j (r,t)*he(v,s)} = JJ e ~ 27ri ^' 1T ~ Xv ^'K {rjj (i, 7)*J#(s, A)} d 7 dA = 0, 
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and, when j = £, 

E{h j {T,t)*h j {v,s)}= II e- 2 ^ T - x ^C J {t, 1 )5{t-s)8{ 1 -\)d 1 d\ 



= P hj ( T -v,t)S(t-s), 
where -P/^r, t) = J e^ 2n ' nT Cj(t,^) d 7 . Moreover, as a result of ([5]), we have 

Ph{T,t)= [ e - 2 ^C(t, 7 )d 7 



J e - 2 ^Cj(t - flj T,7 - b j B)d 1 



./ 

= ^2e- 2 ^ B P h] (r,t-a 3 T), 
i=i 

Next, we derive a formula similar to @ for the scattering functions Cj(t,j). 
Lemma 2. With rjj, hj, P^. and Cj defined as above, we have 



where r is arbitrary. 

Proof. By equations (HJ and ([6]), we have 



E{»fc(*,7)*»&(»,A)} = i E e 2 ^ + ^-^) E {^(r+f,0*^(|, S )} 

n,m£Z 

= ^ E e-<^-^( T+ Ii^,t)*( t - a ) 

1 n 



1 n 



(0) 



(7) 



B 

Since E{77j(t,7)*^ (s, A)} = S(t - s)5( 7 - X)Cj(t, 7), the result follows. □ 

We are now ready to prove the main result of this paper. 

Proof of \Theorem 1[ Let cj, be a J-periodic complex sequence. Transmit x{t) — ^2 ke x Ck ^ kT (^) anc ^ 
observe echo y(t) = Hx(t). Also, define 

y n {t) = (fTs)(i + »T) = H (J2 c khr) (t + nT) = Y^ c k h(t + nT, t + (n - k)T), (8) 

k k£Z 

for < t < T and n 6 Z. Clearly, if h(t,r) is a stochastic process, so are y„(f). In fact, for 
< t, t < T and n € Z, we have 



E {j/n(i)*y(r)} = E C ^ E W* + nT ' * + ( n - k ) T )* h (r, t - £T)} 
fe.zez 

= 4c e P h {nT + t-r,t + (n- k)T)6(t -r-(n-k + i)T). 
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But for < t, t < T, 

S(t-T -{n-k + t)T) 

Thus 



6(t-r), tfn-k + £ = 0, 
0, if n - k + t + 0. 



E {y n (t)*y(T)} = c* k c k - n P h (nT, t + (n - k)T)S(t - r) 
feez 

= Y, c l +n CkPh(nT,t-kT)S(t-T) (9) 

fcez 

= Il n (t)S(t-T), 

where we defined 

IL n (t) = J24+nCkPh(nT,t-kT). 

feez 

Using (JTJ) and ([9]), we obtain 

Ii r+m j(t) = 4+r+mjCkPh (§ + rT, t - fcT) 

feez 

= E «S+r** E (g + rT,t - (fc + a,)T) 

fcez ?=i 

= E e- 2 ^ /J E c *+r c *^i (f + rT, i - (fc + 0i )T) . 
j=i fcez 

Recall Ji3T = 1 and denote for r = 1, . . . , J 

Sr(*,7) = \ E e 2 ^ T ( mJ+ ^n mJ+r (i). 

By Lemma [2] we compute 

fiv(*,7) = \ E e 2 ^ T( " lJ+I ' ) n mJ+r (t) 

m6Z 

= E e - 2 " rb3/J E c *-H-<*£ £ e^^W)^. ((mJ + r)T, t - (fc + a,)T) (10) 

j=l fcGZ mez 

,7 

= E e" 2 ^/ J E 4 +r c fe Q(t - (fc + a,)T, 7 ). 
i=i feez 

But, Cj(t, 7) ^ only if < t < T. Then, C 3 (t - kT - a^T, 7) ^ if and only if 

i-T< (fc + aj)T<t. 

Since < t < T, we must have that Cj(t - (fc + a,)T, 7) ^ only if fc + Oj = 0. Therefore, by ([TO]). 

S r (t, 7 ) = E e- 27ri ^ /J <_ .c_ 0j C j (t, 7 ). (11) 
j'=i 



Let {ck,i} J k i=Q be the Weyl-Heisenberg frame for <C J , where 

c k ,i(r)= e~ 2 ™ l / J c^ k , r = l,...,J. 

The coefficient sequence {ck] k=0 can be chosen so that any J element subset of the so-called Gabor 
frame vectors is linearly independent [KPR08 , a condition of the frame known as the Haar property. 
Moreover, the selection procedure allows us to choose c to be unimodular. In fact, choosing the 
entries of c randomly from a uniform distribution on a unit circle guarantees the Haar property will 
hold with probability one. In particular, this means that a matrix of our interest, determined by 
the geometry of the support set of C(t,j), is invertible. Let's denote by U the J x J matrix whose 
j-th column is Ca jt bj, that is, 

U — \pai ,bi ^025^2 I ' | ^ctj ,bj 1 ■ 

Then, (111 p can be expressed in the matrix form as 

[5 r (t,7)]^ 1 = Cr*[ c _ 0i Cr J (t, 7 )]/ =1 . 
so that we can rewrite the above equation as 

[s r (tn)] J r=1 = u*D[c :i (t n )}^ 1 

with D G C' ,xj a non-singular diagonal matrix with elements of c on the diagonal, that is, D^j = 
c^ aj 5ij. Let V — [Vj r ]' 7 J -_j be the inverse of U*D C . We can now recover Cj(£,7) pointwise from the 
autocorrelation of the output signal with an explicit reconstruction formula derived from (|5]), © 
and Unj), that is, for any j e [1, J] and t € [0, T), 7 e [0, B) 

J J 

C;(*.7) =J2V jr S r (t,7) =Y, V *b ^ e 2 ^ T (™ J+ '')n mJ+r (i). (12) 

r=l r=l meZ 

where II n (t)5(£ - r) = E + nT)*y(r)} . This completes the proof of lTheoremT] □ 

4 Scattering function estimation 

Since h(t,r) is a stochastic process, so is the returned echo y(t). Although ITheorem 11 enables us to 
reconstruct the scattering function perfectly from the stochastic process y(t), in practice we do not 
know the values of E {j/„(t)*y(r)}. Secondly, the delta trains that we use for identification of the 
channel have infinite energy, hence, so do their echoes. To address the first problem, we sound the 
channel repeatedly with the same input signal x{t) = ^2kez c k^kr(t) to obtain y^'(t), . . . ,y( L '(t), 
L independent identically distributed samples of the channel output. As the channel is WSSUS, 
repeated soundings do not pose an insurmountable problem in general. Since the primary purpose of 
this paper is a proof of concept for the reconstruction of the scattering function with non-rectangular 
support, we do not explore the effects of time-gating the delta train. 

Regarding the second obstacle of infinite energy responses, we consider as data for our estimator 
the functions H n (t) defined in (0) for t, r £ [0, T] as follows. 

U mJ+r (t)S(t-T) = i£yW( t _( ro j + r ) T )V')( T ). 
i=i 

Clearly, plugging the above into (fT2)) this induces an unbiased estimator C(t, 7) of the scattering 
function C(t, 7). 

r=l m£Z 
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We need the following trivial observation. We renumber the delta train 

x(t) = ^ c k S ( f -(nJ -mJ + k-r) 
fee [J] n&Z 

and observe 

y(t- (mJ + r)T) 
= jj V{r, 7) e 27ri ^ mJ+r ^x{t - (mJ + r)T - r) dr d 7 

= c kJ2 if e 2lTt(t - (mJ+r)Th S(t - (mJ + r)T - (nJ - mJ - k - r)T - r) dr d 7 

fe€[J] nGZ 



E c feE / #-M-fc)7 1 ,7)e 27rt7( *- ((mJ+nJ - fc+r)T) d 7 
E » ? (i + jfeT,7)e 2 ^< t -< mJ - fc+r > T > d 7 , 



fee [J] 

where n now depends on fc, m and r, and fc, as well as the sequence c r k , depend on r. 

Let us pre-compute the covariance between H mi j +r (t) and U m . 2 j +r (t). As different measurements 
j/O and ?/' ) are independent, it is enough to consider 

e n mi J+r (h)ii m2 J+r (t 2 )*5(ti - 1 3 )5(*2 - * 4 ) 

= ^E|y(i! - (mi J + r)T)*y(t 3 )y(t 2 - (m 2 J + r)T)y(t 4 )^ 

= \ J2 c ii c k2Ck 3 c* ki 

k 1 ,k 2 ,k 3 ,k 4 ,£[J] n 1 ,n 2 ,n3,n 4 e2 

E{?7 (tx - (m J - fci)T, ^)* r? (t 2 - (n 2 J - k 2 )T, v 2 ) 

x 77 (t 3 - {n 3 J - fc 3 )T, 1/3) r? (t 4 - (n 4 J - fc 4 )T, ^ 4 )* } 

x e -27ri[(ti-(miJ+niJ-fci+r)T)i/i+(t4-(n 4 J-fc4)T)i/ 4 ] 

e27 ri[(t3-(n 3 J-fe3)T)i/3+(t2-(m 2 J+n2J-fe2+r)T)i/ 2 ] ^ ^ ^ ^ 

To estimate the variance, we assume that all r](t, 7) are jointly proper Gaussian processes. With the 
help of the Isserlis' moment theorem, we simplify the expectation 

E{ V fa - (m J - h)T, v x )* r] (t 2 - (n 2 J - k 2 )T, v 2 ) 

V(h ~ (n 3 J - k 3 )T, v 3 )t](t 4 - ( n 4J - k 4 )T, v 4 )*} 
= C(ti - (m J - h)T, Vl )S(h - (m J - h)T - (t 2 - (n 2 J - fc 2 )T))<5(^ - v 2 ) 

x C(t 3 - (n 3 J - k 3 )T, u 3 )8{t 3 - (n 3 J - k 3 )T - (t 4 - (n 4 J - k 4 )T))5{u 3 - v 4 ) 
+ C(ti - (m J - k x )T, u)5(t! - (m J - fci)T - (t 3 - (n 3 J - fc 3 )T))<f(i/i - ^3) 

x C(< 2 - (n 2 J - k 2 )T, v 2 )8(t 2 - (n 2 J - k 2 )T - (t 4 - (n 4 J - fc 4 )T))^ 2 - 1/4) 
+ C(ii - (m J - fci)T, - (m J - fci)T - (i 4 - (n 4 J - fc 4 )T))<5(z/i - v 4 ) 

x C(i 2 - (n 2 J - fc 2 )T, ^ 2 )<5(i 2 - (n 2 J - fc 2 )T - (t 3 - (n 3 J - k 3 )T))8{u 2 - v 3 ). 

All three summands can be estimated similarly. We provide detail for the first summand only. Con- 
siderations of delta functions and support constraints of the scattering function cause cancellations. 
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In particular, we have n± = n 2 = 0, k\ = k 2 , n 3 = n 4 = 0, k 3 = k±, 
h{t)= E c *k x c k^c* ki E jjjj dz/i dz/ 2 di/ 3 dz/ 4 

k 1 ,k 2 ,k 3 ,k 4 e[J] m,n2,ns,m& 

x C(h - (m J - fci)T, - (m J - fci)T - (t 2 - (n 2 J - k 2 )T))5(v L - v 2 ) 

x C(t 3 - (n 3 J - fc 3 )T, i^)«5(t 3 - (n 3 J - h)T - (t 4 - (n 4 J - fc 4 )T))<K*/ 3 - v*) 



x e 



2ni[-(t 1 -(m 1 J+n 1 J-k 1 +r)T)v 1 +(^t2-(m2J+n2J-k2+r)T)v2+(t3-(n 3 J-k3)T)v3-{t4-(n i J-k4)T)v4] 



E c ki c *i c *s c Jls E JJCfa + kiT^e-^-^Cfa + ksT,!*) 
k 1 ,k 3 e[J] m,nsez 

-2«(m 1 -m ! )JT, 1 ^ ^ 



x e 



Hence, we have for the covariance 

cov (n miJ+r (t),n m2 , 7+r (t)) - ii(t) + / 2 (t) + J 3 (t). 

We can now estimate the variance on (t, 7) G [0, T)[0, 5). 



|var S r (t, 7)| 
1 

~ B2 



1 

B2 



E e 2m7T J(mi - m2) cov n mi , 7+r (*)n ma , 7+r (t) : 

u,m2 

E e 2 ^ TJ < mi - m2 ) (Ji(t) + I 2 (i) + J 3 (t)) 

Hi">2 

^2 I E l^l 2 / C(t + fciT,i^)X; e™«-^ TJm > A Vl 

feiG[J] mi 
fe 3 G[J] ™ 2 

E |c fcl I 2 / C(t + fciT, ^) E e 2 ^ +v ^ JT+rT ^ dvt 
2 e[./] " ">2 

E ( C ^i) 2 / C(t + fclT,^)E e 2 -(7-2-i)TJm le -4^ l( t +fel T-rT) dz/i 
E 4 / C(t + k 2 T, V 2 ) E e -2-(7-^)T./m 2e 4™ 2(t -fe 2 T- r T) ^ 



feie[J] 
x 



fei6[J] 

X 



fc 2 e[J] 
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1 

TP 



( ]T \c k \ 2 C(t + fcT, 7 )) + ( ]T \c k \ 2 C(t + kT, B - 7 ) 
fce[J] fee[J] 

2 (c/c) 2 (C(t + kT, 7 /2) + C(t + kT, (B + 7 )/2)) ' 2 
fee [J] 



^ E + fcr, T )| 2 + + ^r, i? - T )| 2 

fee [J] 

+ 2 |C(i + fcT, 7 /2)| 2 + 2 |C(t + fcT(S + 7 )/2)| ; 



This gives us a pointwise estimate. We now average to obtain the following estimate for the 
variance of Cj, 



var Cj 




btJJ 







dt d 7 < 



LB 2 



This completes the bias- variance analysis. 



5 Conclusion 

In this paper, we have shown that it is possible to recover the scattering function of a target using 
the second order statistics of the returned echoes from sounding by a custom weighted delta train, 
provided that supp C{t, 7 ) is contained in a union of rectangles of total area 1. This includes channels 
commonly considered overspread under the classic criterion that the area of the bounding box is less 
than 1. In the abstract setting of a continuous WSSUS channel sounded by an infinite delta train, 
being given complete second-order statistics of the received signal, we guarantee an exact recovery 
of the scattering function. It is interesting to note that a single signal suffices for recovery of an 
entire collection of channels which share the gridding of the support set of their scattering functions, 
but not the support itself. We suggest that the proposed technique to reshuffle the patches of the 
support of the scattering function can and should be seen as a replacement of the Zak transform in 
the efficient implementation of the averaged periodogram estimator procedure offered in [AMH99 . 

For completeness, we give an explicit recipe of a simple unbiased estimator, with a predictable 
behavior of the variance given repeated soundings under the assumption of pointwise jointly proper 
Gaussianity of the spreading function. It would be interesting to explore the stability of the technique 
under the time-gating of the input signal. We currently do not pursue this direction, since the 
purpose of this paper is to establish the possibility of the scattering function reconstruction and 
to demonstrate the usefulness of the weighted delta train technique in the field of radar and sonar 
acquisition. 
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